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Abstract 

Two-body relaxation times of nuclear star clusters are short enough that gravitational encounters should 
substantially affect their structure in 10 Gyr or less. In nuclear star clusters without massive black holes, 
dynamical evolution is a competition between core collapse, which causes densities to increase, and heat input 
from the surrounding galaxy, which causes densities to decrease. The maximum extent of a nucleus that can 
resist expansion is derived numerically for a wide range of initial conditions; observed nuclei are shown to 
be compact enough to resist expansion, although there may have been an earlier generation of low-density 
nuclei that were dissolved. An evolutionary model for NGC 205 is presented which suggests that the nucleus 
of this galaxy has already undergone core collapse. Adding a massive black hole to a nucleus inhibits core 
collapse, and nuclear star clusters with black holes always expand, due primarily to heat input from the galaxy 
and secondarily to heating from stellar disruptions. The expansion rate is smaller for larger black holes due to 
the smaller temperature difference between galaxy and nucleus when the black hole is large. The rate of stellar 
tidal disruptions and its variation with time are computed for a variety of initial models. The disruption rate 
generally decreases with time due to the evolving nuclear density, particularly in the faintest galaxies, assuming 
that scaling relations derived for luminous galaxies can be extended to low luminosities. 

Subject headings: galaxies: evolution - galaxies: nuclei - galaxies: dynamics 



1. INTRODUCTION 

Many galaxies have compact stellar nuclei. A nearby 
example is NGC 205, the dwarf elliptical (dE) companion 
to the Andromeda galaxy. NGC 205 exhibits a sharp up- 
turn in its surface brightness profile inside of a few arcsec- 
onds, attributable at l east in part to a population of youn g 
stars near the center (iHodgd 119731: ICarter & Sadler! I1990I) . 
The mass of the nucleus is ~ 10 6 M W and its radius is ~ 
IOp c (iJones et d]ll996tjVglluri et alJl2005h . iReavesi (1 19771) 
and iRomanishin et al.l (1 19771) found similar central light ex- 
cesses in many ear ly-type galaxies in the Virgo cluster; 
Bingg erTet al.1 (11987). in their comprehensive survey of Virgo, 
found that roughly one-fourth of the dwarf (dE and dSO type) 
galaxies were nucleated. The nuclei in the Virgo galaxies are 
sufficiently compact that they appear unresolved in ground- 
based images, and their detection in these photographic sur- 
veys was facilitated by the low central surface brightnesses of 
dwarf galaxies. 

Recent observations with the Hubble Space Telescope have 
revealed that nuclear star clusters are present in a much wider 
variety vof galaxies, including lu minous elliptical galaxies 
and the bulges of spiral galaxies dCarollo et al.lll997l 119981 ; 
iBoker et alJT2002t ICote et al.ll2006l) . The fraction of galaxies 
with clearly distinguishable nuclear clusters is now believed 
to be as high as 50%-75%. The nuclei are not present in E 
galaxies brighter than absolute blue magnitude Mb w — 19, 
but this may be due to the steeply-rising luminosity pro- 
files in these galaxies which would make detection of a nu- 
clear excess difficult. The frequency of nucleation also ap- 
pears to drop in the faintest galaxies, falling essentially to 
zero at galaxy abso l ute magnitudes Mr ss —12 and below 
dSandage et alj|1985t Ivan" den i Berghlll986l) . 

Nuclear star clusters have the following properties dBokerl 
l2007t Ivan der Marel et ai1l2007l and references therein), (a) 



They are 10-100 times brighter than typical Milky Way glob- 
ular clusters, (b) Their sizes correlate with their luminosities 
as r ~ L° 5 . (c) Spectra reveal extended star formation histo- 
ries, and the luminosity-weighted mean stellar age correlates 
with galaxy Hubble type, from ~ 10 Myr in late-type spirals 
to ~ 10 Gyr in ellipticals. However the mass is typically dom- 
inated by an old stellar population, (d) Nuclear masses corre- 
late with host galaxy (bulge) masses as M nuc « 0.003Mb u i ge , 
as well as with other global properties of the bulge such as 
velocity dispersion and Sersic index. 

The fact that nuclear star clusters obey similar scaling 
relations with host galaxy properties as do supermassive 
black holes has led to the suggestion that there may be a 
fundamental connection between the two cla sses of object 
dFerrarese et al.ll2006al;IWehner & Harrisl2006l) : perhaps both 
should be assigned to the single category of "central massive 
object". This suggestion raises the question of whether nu- 
clear star clusters can co-exist with supermassive black holes. 
ISeth et all (120081) list roughly a dozen galaxies with nuclear 
emission lines or X-ray luminosities indicative of an active 
nucleus, and which also contain nuclear star clusters. 

This paper addresses the evolution of nuclear star clusters 
due to random gravitational encounters between stars. Ignor- 
ing the possible presence of a massive black hole, relaxation 
times in the nuclei work out to be about an order of magnitude 
longer than those in globular clusters, but still short enough 
that gravitational encounters would substantially affect their 
structure after 10 Gyr. For instance, in NGC 205, which has 
one of the best-resolved nuclear star clusters, the central re- 
laxation time is only ~ 10 8 yr, short enough that core collapse 
could have already occurred (IJones et al.lll996t IValluri et al.l 
120051) . 

As p oint ed out by Dokucha ev & Ozernoil (1 19851) . iKandrupl 
( 1990) and lOuinianl (1 19961 ). evolution of a star cluster at the 
center of a galaxy can be qualitatively different than that of an 
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isolated stellar system. If the nucleus is sufficiently diffuse, 
its velocity dispersion will be lower than that of the surround- 
ing galaxy, and the usual tendency toward core collapse will 
be opposed by a flow of energy from galaxy to nucleus. The 
fate of the nucleus - collapse versus expansion - depends on 
whether the core collapse time is short or long compared to 
the time for th e nucleus to absorb energy from the galaxy. 
lOuinlanl |[l996) analyzed a simple family of two-component 
models and derived a condition for the minimum compactness 
of a nucleus in order to resist expansion. Quinlan's condi- 
tion implies a maximum size r nuc /r ga i « 0.02 for nuclei with 
masses M nuc /M Ba i = 10~ 3 . This is intriguingly close to the 
observed sizes of nuclear star clusters. 

Even a nucleus compact enough to satisfy Quinlan's crite- 
rion could escape core collapse if it contained a black hole. 
A black hole al so acts as a h eat source, by destroying and 
absorbing stars (Sha pirolll977l) . Nuclear star clusters contain- 
ing massive black holes should therefore expand, regardless 
of their initial degree of compactness. 

In this paper, integrations of the orbit-averaged Fokker- 
Planck equation are carried out for models of nucleated galax- 
ies with and without black holes. The Fokker-Planck model 
includes a loss term describing scattering of stars into the tidal 
disruption sphere of a black hole, if present. Physical colli- 
sions between stars are ignored, and the mass liberated from 
disrupted stars is assumed to escape instantaneously from the 
galaxy; these are both approximations but not unreasonable 
ones (e.g. in NGC 205 the timescale for physical collisions 
is ~ 3 orders of magnitude longer than the relaxation time). 
Initial conditions are derived from mass models that mimic 
the observed luminosity profiles of nucleated galaxies. Quin- 
lan's criterion for black-hole-free nuclei is refined; we find 
that observed nuclei are either in the core-collapse regime, 
or have sufficiently long relaxation times that they would not 
have evolved appreciably in 10 Gyr. In models containing 
black holes, the dominant process driving nuclear evolution is 
almost always found to be heating by the galaxy. 

The models presented here are idealizations in one impor- 
tant respect: only a single stellar mass group is considered 
and stellar masses are assumed to be fixed in time. The ef- 
fects of mass segregation, star formation and stellar evolution 
are therefore ignored. Multi-mass models will be considered 
in a later paper. While the calculations presented here have 
some relevance to nuclear star clusters in all galaxies, the fo- 
cus is on models that are structurally similar to low-luminosity 
elliptical galaxies. 

Section II summarizes the properties of nuclear star clus- 
ters and their host galaxies that are relevant to this study. Sec- 
tions III and IV describe the results of evolutionary calcula- 
tions without, and with, black holes, respectively. Section V 
sums up. 

2. NUCLEAR STAR CLUSTERS AND THEIR HOST GALAXIES 

2.1. Nucleated Galaxies in the Virgo Cluster 
FigureQJ shows half-light radii r e g- and masses M nuc of nu- 



clei o bserved in the ACS Virgo Cluster Survey dCote et alj 
2004). Dashed lines indicate constant values of the half-mass 
relaxation time 

1.7xl0 5 [r A (pc)] 3 / 2 ArV2 

f * = r 2 y ears < 1 > 

[m/M & } 1/2 

( Spitz erl 19871) . with r/, set equal to r e g. The mean stellar mass 
m was set to one solar mass and N was computed from N = 



Mauc/m. Half-mass relaxation times for many of the nuclei 
are less than 10 Gyr. 

Figure QJ> shows the dependence of relaxation time on 
galaxy absolute blue magnitude Mb- Relaxation times are 
plotted both for the nuclei, and for their host galaxies; in the 
case of the galaxies, r/, in equation (Q~|) was equated with the 
observed effective radius R e (e.g. Ciotti 1991) and N was set 
to Mg a i/ra. Dotted lines are least-squares fits to logf r /, vs. Mb 
for galaxies and nuclei independently; the fitted relations are 

logio(%ai/yr) = 14.2 - 0.336(M fi + 16), (2) 
logio(Wyr) = 9.38 - 0.434(M fi + 16) (3) 

where we have written f ga i and f nuc for the half-mass relax- 
ation time of the galaxy and the nucleus respectively. The 
second of these relations predicts nuclear relaxation times be- 
low 0. 1 Gyr in galaxies with Mb ~ —12, the least lumino us 
galaxies that contain compact nuclei ( Ivan den Be rgh 1986). 

Figure [TJ; plots nuclei in the size ratio vs mass ratio plane. 
The symbol sizes in this plot are proportional to log 10 f nuc ; 
filled symbols have t nuc < 3 Gyr and open symbols have 
fnuc > 3 Gyr. There is evidence in this figure for two pop- 
ulations: compact nuclei with short relaxation times and 
more extended nuclei with long relaxation times. Figure d, 
combined with the k nowledge that nuclear masses correlate 
with g alaxy masses dWehner & Harris! [20061 : iFerrarese et alJ 
2006a), impl ies that the large-f nuc population in Figure [TJ; is 
associated with the brightest galaxies. The lack of nuclei with 
r e g/R e <, 10~ 3 probably reflects an observational limitation: 
at the distance of Virgo, the ACS/WFC pixel scale of ~ 0.049 
arc second corresponds to ~ 3.9 pc, or ~ 0.004 times a typical 
galaxy effective radius of ~ 1 kpc. 

2.2. NGC 205 

Nuclear star clusters beyond the Local Group are too dis- 
tant for their structure to be well resolved. The nearby (~ 824 
kpc; McConnachie et al. 2005) dE galaxy NGC 205 con- 
tains one of the best resolved nuclei. Figure [2] shows surface 
brightness data of NGC 205 in the I band, obtained by com- 
bining HST ACS ob servations fromlValluri et alJ d2005l) with 
ground-based data of lKim & Leeldl998l) and lLeeldl996l) . The 

1 /2 

horizontal axis is tsma [1 — E(r)J ' where tsma is the isopho- 
tal semi-major axis and e is the ellipticity; NGC 205 is mod- 
erately elongated, (e) « 0.5. 

These circularly-symmetrized data were fit to a two- 
component spherical model defined, in terms of the (3d) lu- 
minosity density j(r), as 



Kr)=Mr)+j 2 (r) 
j2(r)=j n uc(l+r 2 /r 2 c y y/2 . 



(4) 

(5) 
(6) 



The galaxy is represented by Einasto's (1965) law, identical 
in functional form to the law of Sersic (1963) that is more typ- 
ically applied to surface density profiles. The quantity bin) 
is chosen such that r\n is the (3d) radius containing 1/2 of 
the total mass corresponding to j\\ n is an index that con- 
trols the curvature of the profile (not to be confused with the 
Sersic index; the latter is typically smaller by ~ 1 — 2 than the 
Einasto n that describes the space density). The nucleus is rep- 
resented by a generalized "Hubble" law; r c is the core radius 
and y controls the steepness of the nuclear density falloff, e.g. 
y = 5 is Plummer's (1911) model and y = 2 is the isothermal 
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FIG . 1 . — Prop erties of nuclear star clusters in the subset of Virgo cluster early-type galaxies that were found to be nucleated in the ACS Virgo Cluster Survey 
I Cote et al. 2004). (a) Nuclear radii and masses; masses are from the compilation of Seth et al. (2008). Dashed lines correspond to nuclear half-mass relaxation 
times of (10 s , 10 9 , 10 10 , 10" , 10 12 ) years increasing up and to the right, (b) Half-mass relaxation times of nuclei (stars) and their host galaxies (filled circles). 
Dotted lines are least-squares fits, (c) Vertical axis is the ratio of half-light radii of nuclei (r c ff) to half-light radii of galaxies (R e ); horizontal axis is the ratio of 
nuclear mass to galaxy mass. Symbol size is proportional to the logarithm of the nuclear relaxation time in (b). Open circles have logio(f r /,/yr) > 9.5 and filled 
circles have log 10 (f r ft/yr) < 9.5. Dashed lines indicate the critical value of r e g/R e above which nuclei expand, for galaxies described by Einasto indices n = 2 
(top), n = 3 and n = 4 (bottom). 

two-component model of equation ( 12.2b will be used as a ba- 
sis for the evolutionary studies described below. We note 
that one parameter in this fit, the Einasto index n (or rather, 
the corresponding Sersic index describing the surface bright- 
ness profile; Einasto fits have not yet been carried out for 
many galaxies) is known to correlate roughly with galaxy 
luminosity, in the sense t hat brighter galaxies have larger n 
dGraham & Guzmanll2003l) . 
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FIG. 2. — Two-component (eq. 12.21 lit to surface brightness data of 
NGC205. Thin solid lines show the nuclear and galactic components sep- 
arately, after convolution with the ACS point-spread function; dashed line is 
before PSF convolution. Lower panel show surface-brightness residuals, n: 
Einasto index of the galaxy component; r c : nuclear core radius; y: power-law 
index of the nuclear space density. 

sphere. Equation (I2.21 i was compared to the data of Figure [2] 
after spatial projection and convolution with the circularly- 
symmetrized ACS HRC point-spread function. The optimiza- 
tion routine adjusted the five parameters (jgai,jmici r i/2, r c,n) 
at fixed y; the latter parameter was varied in steps of 0.1. 

Figure|2]shows the fit that minimizes the surface-brightness 
residuals, with y = 3.2. The fit is good (~ 0.05 magnitude 
rms deviation), although this decomposition of the brightness 
profile into "galaxy" and "nucleus" is not certainly not unique. 

Since it provides such a good description of NGC 205, the 



3. NUCLEAR EVOLUTION WITHOUT BLACK HOLES 

We first consider evolution, due to gravitational encounters, 
of a galaxy containing a compact nucleus consisting of stars 
with a single mass and no nuclear black hole. Physical colli- 
sions between stars are ignored. In this idealized model, evo- 
lution takes place on a time scale determined roughly by the 
central (nuclear) relaxation time, t rc ; stars are scattered into 
less-bound orbits, transferring heat to the nuclear envelope 
and causing the central density to increase. If the nucleus were 
isolated, and ignoring sources of heat like hard binar y stars, 
core c ollapse would take place in a time of ~ 10 2 f re ( Spitzer 
1 1987b . However because the nucleus is embedded in a galaxy, 
flow of heat can also take place in the opposite sense, from 
galaxy to nucleus, causing the latter to expand; core collapse 
is the n delayed until a much longer time has elapsed (Ouinlan 
1996). 



3.1. Evolution timescales 

We begin by deriving an approximate expression for the 
time scale to transfer energy from a galaxy to a nucleus, 
then derive an approximate condition for heating to reverse 
core collapse. Our approa ch i s similar to tre atments in 
iDokuchaev & Ozernoil (11985b and lKandrupl (Il990h . 

Consider a two-component system consisting of galaxy and 
nucleus. Both components are assumed homogeneous and are 
characterized by a density (p ga i,pnucX radius (r ga i,r nuc ), rms 
velocity (V ga i, Vnuc), and half-mass relaxation time (f ga i,f nU c)- 

Define e nuc = (l/2)p nuc V n 2 uc to be the kinetic energy per 
unit volume of the nucleus. Assuming Maxwellian velocity 
distributions and a single stellar mass m, the rate of change of 
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FIG. 3. — Density (upper) and velocity dispersion ( lowe r) profiles of a set 
of galaxy models constructed according to equation )2.2t . normalized such 
that G = M ga i = 1. All models have M llllc /M g!1 i = 0.003, n = 3 and y — 
4. Different curves correspond to different degrees of nuclear compactness, 
i.e. logio(r nuc /r ga i) = (-3.72,-3.20,-2.72,-2.20,-2.02,-1.90,-1.72). 
Thick lines denote models in which the nucleus undergoes "prompt" core 
collapse; thin lines are models in which heat transfer from the galaxy causes 
the nucleus to expand initially (as show in Fig. [4). 



e nuc due to gravitational encounters is 



ds n 



dt 



4V67T.G mpnucPgallnA 



(V 2 ,-V 2 "l 

l * gal * nuc I 



-v 2 



(7) 



(Spitzer 1987, eq. 2-60). A necessary and sufficient condition 
for the nucleus to be heated by the galaxy is V ga \ > V nuc . When 
this condition is satisfied, the nuclear heating time is 



^heat — 



1 c/e n 



dt 



-l 



(8) 



Pnuc 

48 V JtG 2 mp nuc lnA V Pgai 



„ \ 3 / 2 

(v g yv n 2 uc +i 

(^gal/^nuc 



-(9) 



In the limiting case V„^i ^> V n uc this becomes 



the 




(10) 



77^ Ual (ID 



(12) 



Models that satisfy the condition V ga i > V nuc have 
Pnuc <^ Pgai- Thus, the nuclear heating time is of the same 
order as, or somewhat less than, the geometric mean of f nuc 
and fg a i. According to Figure lb, this time is shorter than 10 
Gyr in at least some galaxies. 

Heating from the galaxy will reverse core collapse if T neat is 
shorter than the nuclear core collapse time. One definition of 
the latter time is 



1 dp n 



dt 



-l 



(13) 



where varies from ~ 10 in the early stages of core col- 
lapse to an asymptotic value of ~ 300 (e.g. Cohn 1980). The 
condition T ne at < ^cc becomes 



Vgal 



Pgai V m 



-1 



(14) 



where V ga i ^ V nuc has agai n been assumed. This expres sion is 
similar to equation (14) of iDokuchaev & Ozernoil d!985l) . We 
can convert equation (fl4l > into a relation between the quanti- 
ties plotted in Figure 1 by writing p nU c ~ M nuc /rl uc and by 
applying the virial theorem separately to both components, 
i.e. 



V 1 



GM n 



V, 



GM. 



gal 



gal _ 



(15) 



the former expression will only be approximately true for low- 
density nuclei. With these substitutions, the condition ( fl4] > 
becomes 

:, 10 ; f'i-* ' 
U00, 



Mi 



gal 



'gal 



(16) 



Setting Mnuc/^gai = 0.003, equation ( fT6b implies that 
rnucAgai must be smaller than ~ 0.05 in order for core col- 
lapse to occur, a value that lies within the range of observed 
nuclear sizes (Fig. lc). While the numerical coefficients in 
equation ( fl"6l ) should not be considered accurate, the general 
form of the relation will be verified below via more realistic 
evolutionary models. 

3.2. Fokker-Planck model 

Evolution was modelle d using the isotropic, orbit-averaged 
Fokker-Planc k equation (jCohnl 19801) . in the implementation 
described by iMerritt et alj (1200*71) ! The initial galaxy was 
given a density law as defined by equation ( 12.2b ; the self- 
consistent gravitational potential and isotropic phase-space 
distribution functions were derived from Poisson's and Ed- 
dington's formulae respectively. Henceforth we characterize 
these initial models by the two dimensionless parameters 



>"nucAgal, M nuc /M gai 



(17) 



where r nuc and r ga i are the spatial (not projected) half-mass 
radii of nucleus and galaxy respectively, and M nuc and M„ a \ 
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FIG. 4. — Evolution toward core collapse in the galaxy models plotted in Fig. [5] Vertical axis is the core density normalized to its initial value. The horizontal 
axis is the time in units of the initial, half-mass relaxation time f,;,, defined by the galaxy (left) or by the nucleus (right). 



are the total masses corresponding to the two components. 
Given the varied and non-unique ways in which half-light 
radii are defined observationally, we will henceforth simply 
identify r n ucA ga i with r e f[/R e when comparing models to data. 
(This identification should probably only be trusted to within 
a factor of two or so.) 

Figure [3] shows p(r) and a(r) for a set of galaxy mod- 
els with Einasto index n = 3 and nuclear slope y = 4; a(r) 
is the Id velocity dispersion. Each of these models has 
Mnuc/Mgai = 0.003 but they differ in their degree of nuclear 
concentration, from r nU c/r ga i ~ 0.0002 to ~ 0.03. The de- 
pendence of o on r is complex, sometimes exhibiting mul- 
tiple maxima; but as the nucleus is made less compact, the 
central velocity dispersion drops, and when r nuc /r ga i <; 0.003 
the central a is lower than the peak value in the outer galaxy. 
"Temperature inversions" like these imply a flow of heat from 
galaxy to nucleus, which tends to counteract the outward flow 
that would otherwise drive th e nucleus toward core collapse 
dDokuchaev & Ozernoilll985l) . 

Figure|4]shows the change with time of the central densities 
of the models in Figure [3] Times are normalized by the ini- 
tial, half-mass relaxation time; the left panel normalizes to the 
galactic relaxation time f ga i, and the right panel to f nuc . There 
is a critical value of r nuc /V ga i above which evolution toward 
core collapse is halted and the nucleus expands. Core col- 
lapse still occurs in these models but on a much longer time 
scale, roughly ten times the galaxy half-mass relaxation time 
(and far longer than galaxy lifetimes ; see F ig. QJ)). The de- 
lay is due, as first shown by Quinlan (1996), to the galaxy's 
need to first reverse the temperature gradient near the center 
by creating a large, flat core. In the models with nuclei that 
contract, core collapse occurs in 15 — 20 times the initial nu- 
clear half-mass relaxation time. 

Experiments like the ones illustrated in Fig. [4] were carried 
out for initial models with a range of r nuc /r ga i and M nuc /M ga i 
values, and with Einasto indices n = (3,4,5). The critical 
size ratio that separates collapsing-core and expanding-core 
models is plotted in Figure lc for the three different values of 



n. The critical values are well fit by relations of the form 



M %a \ V r gal / 

where 

n = 2: A = 300, 5 = 2.70 

n = 3:A = 390, B = 2.45 

n = 4: A = 355, B = 2.20. (19) 

Quinlan (1996) derived a similar criterion for double- 
Plummer-law galaxy models; his relation is most similar to 
the one found here for n = 2, which is reasonable given the 
low central concentration of Plummer models. While these 
results are all based on nuclei with p ~ r~ 4 envelopes, we 
note that half-mass radii are nearly invariant to the structural 
changes that occur prior to core collapse. Hence the crite- 
ria of equation ( fT9l should apply approximately to a much 
broader class of models, including models in which the nu- 
cleus evolved toward core collapse starting from the initial 
conditions adopted here. 

Interestingly, Figure QJ shows that observed nuclei almost 
all lie in the "prompt" core collapse regime; the only clear 
exceptions are nuclei with such long relaxation times ( > 10 
Gyr) that very little evolution would have occurred since their 
formation. 

It would be nice to go one step further and argue that the 
paucity of nuclei above the critical lines in FigureQJ is due to 
dynamical evolution, i.e. that there once existed a population 
of low-density nuclei that underwent expansion and that now 
have such low densities that they are no longer observable as 
distinct components. This argument is valid only if dissolu- 
tion time scales for the nuclei in this region of the plot are 
shorter than ~ 10 Gyr. This is true for some, but not all, pos- 
sible initial conditions. In models near but above the critical 
line, the time for the central density to drop by a factor ~ 2 is 
found to be ~ 10~ 4 f ga i forM nuc /M ga i <J3x 10~ 4 but increases 
to ~ 10~ 3 fgai for the more typical value M nuc /M ga i <, 3 x 10~ 4 
(e.g. Fig.|4|l. The latter time is longer than 10 Gyr in almost 
all galaxies. 
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Finally, we consider the special case of NGC 205 (Fig. [J)- 
On which side of the divide does it lie? Dynamical mod- 
elling of NGC 205 suggests that the nuclear mass-to-light 
ratio is ~ a few times smaller than that of the g alaxy as a 
whole dCarter&Sadleri[l990l IValluri et al.ll200l . To con- 
struct a mass model, we therefore simply decreased the am- 
plitude of the nuclear component in the luminosity profile fit 
of Figure |2] by a factor of 0.3 at all radii. Figure [5] shows 
the projected "surface brightness" of this mass model, i.e., 
— 2.51og£, where T.(R) is the mass surface density at pro- 
jected radius R. The r nuc /V ga i and M nuc /M ga i values for this 
mass model place NGC 205 close enough to the critical lines 
in Figure[TJ: that its fate is difficult to judge. 

Accordingly, a variety of initial models (all with n = 3.04 
and y = 4) were constructed and integrated forward for ~ 10 
Gyr to see if any could match the reconstructed mass profile 
of NGC 205. Figure [5] shows the results of one such exercise: 
this initial model undergoes "prompt" core collapse, achiev- 
ing singular density in a time of ~ 1.0 x 10~ 4 galaxy half- 
mass relaxation times, i.e. ~ 10 10 yr. The final, projected 
mass profile is essentially identical to that of NGC 205 except 
at R <, 0.1" where the data suggest a fiat core; this difference 
is probably not significant however given that the luminosity 
profile is not well constrained at radii smaller than the PSF 
(Alternatively, NGC 205 might not have reached full core col- 
lapse; or some process like binary heating might have caused 
the core to re-expand.) While not conclusive, this and similar 
experiments suggest that the luminosity profile of NGC 205 
is consistent with this galaxy's nucleus having undergone core 
collapse at some time in the last 10 Gyr. 

It is interesting to note that even slight reductions in the cen- 
tral density of the initial model in Figure|5](larger r nuc , smaller 
Mmc) led t° expansion rather than collapse. In other words: 
if its formation history had been only slightly different, NGC 
205 might no longer contain a compact nucleus. 

4. NUCLEAR EVOLUTION WITH BLACK HOLES 

4.1. Models and Assumptions 

Next we consider the evolution of nucleated galaxies con- 
taining massive black holes; the latter are idealized as New- 
tonian point masses with fixed mass M^. Evolution of f(E) 
describing the stars is again followed via the orbit-averaged 
Fokker-Planck equation with evolving potential, including the 
(fixed) potential from the black hole, and with a modifica- 
tion that accounts for the loss of stars th at are scattered into 
the black hole's sphere of tidal disruption dMerritt et alJ2007h 
The latter is assigned a radius of r t where 



n 



vr 



m I 



(20) 



with a stellar radius and r| w 0.844, the value appropriate to 
an n = 3 polytrope. If E = — v 2 /s + i|/(r) is the binding energy 
per unit mass of a star in the combined potential of the galaxy 
and black hole, then the tidal disruption sphere defines a "loss 
cone" of orbits with angular momenta J such that J < Ji c (E), 
where 



Jf c (E) = 2rf Mr t )-E] « 2GM bh r t . 



(21) 



For black hole masses of interest here (My, < 10 M ), r t is 
always greater than the radius of the event horizon, i.e. stars 
are disrupted not swallowed. 
The Fokker-Planck equation now has the form 



dN dF E , . 

a7 = "a# 



(22) 




R/r, 



/2 



FIG. 5. — Open circles show the surface "mass density" of a model of NGC 
205 derived from Fig. f2] the normalization of the nuclear component was 
reduced by a factor of . 3 relative to that of the galaxy to account for the lower 
mass-to-light ratio of the nucleus. Heavy solid line shows a mass model that 
was integrated forward via the Fokker-Planck equation; progressively thinner 
lines show the surface density profile at (0.8,0.98,0.997,0.9995) times the 
core collapse time of ~ 10 Gyr. 

Here, N(E,t) is the distribution of stellar energies; J {E,t) 
is the flux of stars per unit of energy that are scattered in the 
angular-momentum direction into the tidal disruption sphere; 
and FE(E,t) is the e nergy-directed flux. Each of these terms 
is defined in detail in lMerritt et al.1 (120071) . 

Three additional, dimensionless parameters are associated 
with the black hole models: r ( /r ga i, Mbh/M ga i, and N = 
M ga i/m. The first of these can be written 



-^=2.0xKT 9 

'•gal 



r g al 
1 kpc 



10 6 



1/3 



Re 



it- • (23) 



In the "empty loss cone" limi t, the tidal disruption rate varies 
only logarithmically with r, (Light man & Shapiro! 19771) . The 
dependence on r t is steeper if the loss cone is partially full 
(as is the case in these models); nevertheless, here we sim- 
ply fix = Rq and m = M . Furthermore, half-light radii 
for dE galaxies average ~ 1 kpc with a very weak depen- 
den ce on galaxy luminosity (al though with considerable scat- 
ter; iGraham & Guzmara d2003l) ) and so we set 



^=4.5xl0- 9 

'"gal 



10 



10 



n y /3 /io 3 M b h 



gal 



1/3 



(24) 



Equation d24T > gives the first of the three new parameters in 
terms of the other two. 

We relate N in this expression to observable galaxy proper- 
ties by setting m w M Q and expressing M ga i in terms of galaxy 
absolute magnitude. A typical B-band mass-to-light ratio for 
dE galaxies is 2 <, (M/L)b ^ 4 in Solar units (e.g. Seth et al. 
2008) so the relation between Mb and is roughly 



logi Af« 10 -0.4 (M B +18). 



(25) 



In the evolutionary models presented above, was a trivial 
parameter, affecting only how the time unit of the calculation 
(essentially the relaxation time) was converted into years (e.g. 
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eq. [TJ. In models containing a black hole, N also fixes the 
ratio between orbital periods P(E) and the (orbit-averaged) 
time scales for diffusional loss-cone refilling, or 



q(E)=. 



1 



(26) 



Rlc{E) 

ttP(E)R lc {E)- l t-\ (27) 

Here R = J 2 /J C (E) 2 is a dimensionless angular momentum 
variable, < R < 1, with J C (E) the angular momentum of a 
circular orbit of energy E, and ((AR) ) is the diffusion coeffi- 
cient associated with R. For a galaxy with given structural pa- 
rameters (M ga i,r ga i), q °< N~ l . Small/large q (i.e. large/small 
AQ correspond to empty/ful l loss cone regimes respectively 
dLightman & Shapiro! 1 1977h . We adopt the prescription of 
ICohn & Kulsrud] (1 19781) for relating the loss cone flux J to 
the phase-space density / given q. The integrated loss rate, 
/ J (E)dE, is denoted by M. The contributions of the masses 
of tidally disrupted stars either to the galactic potential or to 
the mass of the black hole are ignored. 

Masses of (putative) black holes in dE galaxies are poorly 
constrained, and so several values for the last of the three 
parameters M\,h/M oa \ were tried, in the range 10~ 4 < 
M bh /M gal < 0.03. 

In order to accomodate a black hole, a a slight modifica- 
tion to the equilibrium galaxy models of equation (12.2b was 
required. In an isotropic stellar system containing a central 
point mass, the central density must increase at least as fast as 
r -1 / 2 in order for f(E) to remain nonnegative. The flat nu- 
clear density profile of equation i2.2\ was therefore changed 
toap« r -° 5 dependence inside r w r c . 

In what follows, we fixed the following parameter values: 

y=4, n = 3. (28) 

Setting n = 3 is reasonable given the weak observed de- 
pendence of galaxy Sers ic index on Mb for dE galaxies 
dGraham & GuzmInS 2003). but we note again that our models 
should not be scaled to giant E galaxies which have steeper 
central density profiles. We then carried out four series of 
integrations, starting from initial models with the following 
parameters: 

• Series I: M nuc /M gal =3x 10~ 3 , M bh /M gal = 1 x 10~ 3 . 
The parameter r nuc /r ga i was varied from 2 x 10~ 4 to 
3 x 10~ 2 . The two mass ratios that define this se- 
ries of models are close to the average observe d values 
dFerrarese & Fordll2005HFerrarese et alj|2006ah . 

• Series II: r nuc / r g a i = 3 x 10 4 The nuclear mass was 
varied over the range 1 x 10~ 4 < M nuc /M ga i < 3 x 10~ 2 
and the black hole mass was set to M nuc /3. 

• Series III: M nuc /M ga i = 3 x 10~ 3 , r nuc /r ga i = 3 x 10~ 4 . 
The black hole mass was varied over the range 10~ 4 < 
M bh /M ga i < 3 x 10 " 



-2 



The number of stars in Series I-III was fixed at N = 10 10 , 
appropriate for a Mb rj — 18 galaxy (eq.l25ll. As discussed be- 
low, the choice of N has very little effect on the evolution of 
the density profile, and the results for p(r,f) can robustly be 
scaled to galaxies of different luminosities after adjusting the 
value of fgai. The assumed value of Af does, however, substan- 
tially affect the dynamics of the loss cone and the rate of tidal 
disruptions, and to investigate this dependence a fourth series 
of integrations were carried out: 



-1 1 I I I I III 1 1 I I Mill 1 1 I I I I III 1 1 I I I I III 




r Agal 

FIG. 6. — Evolution of density (upper) and velocity dispersion (lower) pro- 
files in a model from Series I with M nuc /Mgai = 0.003, r nuc /r gli \ = 0.0002 and 
Mbh/Mg a i = 0.001. Thick curves are the initial models and blue curves show 
the final time step; in units of ? eii [, the displayed times are (0,3 x 10~ 6 ,1 x 
10~ 5 ' v in-4 i v in 3 



5 ,1 x 10" 4 ,1 x 10~ J , 0.01, 0.1). 

• Series IV: M nuc /M gal = 3x 10~ 3 

->-3 



>"nucAgal = 1 X 10 3 , 

Mbh/M ga i = 1 x 10~ 3 . The parameter Af was varied 
from 10 6 to 10 ; according to equation (l25l >. this 
corresponds to a range in galaxy absolute magnitude 
-8;>M B ;> -20.5. 

Each integration was continued until a time of roughly 0. lf ga i. 
Using equations ( |25] | and @, this corresponds to ~ 10 14 yr 
for M B = -18, 10 13 yrforM B = -15, 10' 2 yrforM B = -12, 
and 10 11 yr for M B = -9. 

In a galaxy containing a massive black hole, the nuclear 
half-mass relaxation time is not a well defined quantity. In- 
stead, we will cite below the local relaxation time f r (r) in 
these models, defined as 



t r {r) 



0.34o(r) 3 



(29) 



p(r)mG 2 In A 

( Spitz erl [l987l) . where a(r) is the Id stellar velocity disper- 
sion defined above and includes the contribution to the stellar 
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FIG . 7 . — Comparison of the density evolution in two integrations including 
(solid lines) and without (dashed lines) the loss term f representing tidally 
disrupted stars. Initial model (heavy line) hadM nllc /M ga i = 0.003, r nuc /r ga i = 
0.0003 and Mhh/M ga \ = 0.001. In units of f ga i, the displayed times are (0, 1 x 
10~ 5 ,1 x 10~ 4 ,1 x 10" 3 , 0.01, 0.1). 

motions from the presence of the black hole itself. A natu- 
ral radius at which to evaluate t r is the black hole's influence 
radius, 

GM bh 

r D h= — 5~. (30) 
°o 

In this formula, we set CTq equal to its peak value in the black- 
hole-free model; e.g. for the models of Figure[3](n = 3), Go ss 

0. 25. or 

in units where G — M ga i — 1 . This is reas onably equivalent to 
the way that ha is defined by observers. iMerritt et al.l (120071) 
plot f r (H3h) as a function of Oo in a sample of (mostly lumi- 
nous) galaxies; they find t r (r^) m 10 11 yr(ao/100km s -1 ) 7 5 . 

As defined here, rbh is essentially time-independent in these 
simulations, although the values of p(rbh), o(/bh) and hence 
f r('"bh) do evolve substantially. 

4.2. Results: Density Evolution 

Models with central black holes evolve differently than 
models without black holes, in two important ways. 

1. On time scales of order the nuclear relaxtion time, the 
density within the black hole's influence sphere attains the 
p ~ r~ 7/4 "zero-flux" Bahca ll-Wolf quasi-steady-state form 
dBahcall & WolflfT976lfT977h . 

2. The presence of the black hole counteracts both the ten- 
dency toward core collapse in compact nuclei, and the ten- 
dency toward core expansion in low-density nuclei. The black 
hole's mass - while a small fraction of the total galaxy mass 
- is always comparable with the mass of the nuclear cluster in 
these models, and the fixed gravitational potential associated 
with the black hole acts to decouple changes in the central 




10 4 10 3 0.0 



1 gal ' gal 

FIG. 8. — Density profiles at t = (left) and after an elapsed time of 
10~ 4 fgai (right) for various models containing black holes. Line widths corre- 
spond between left and right panels. Top: Series I; r nuc /r ga i = (2 x 10~ 4 ,3 x 
10~ 4 ,1 x 10~ 3 ,3 x 10~ 3 ,1 x 10~ 2 ,3 x 10~ 2 ); line width is a decreasing 
function of r nuc . Middle: Series II; M nllc /M ga i = (1 x 10~ 4 ,3 x 10~ 4 , 1 x 
10~ 3 ,3 x 10~ 3 ,1 x 10~ 2 ,3 x 10~ 2 ); line width is a decreasing function of 
M Imc . Bottom: Series III; M bh /M ga i = 1 x 10~ 4 ,3 x 10~ 4 ,1 x 10~ 3 ,3 x 
10~ 3 , 1 x 10~ 2 ,3 x 10~ 2 ; line width is a decreasing function of Mbh- 

density from changes in the central velocity dispersion, break- 
ing the feedback loop that drives core collapse. Furthermore, 
in low-density nuclei, the presence of the black hole raises the 
"temperature" of the nucleus and decreases the inward heat 
flux from the surrounding galaxy, reducing the expansion rate 
compared with that of a black-hole-free nucleus. 

Because core collapse is suppressed in these models, the long- 
term evolution (f ^ f r (Vbh)) is always toward expansion. 

Figure [6] illustrates the evolution of a model from Series I, 
with M nuc /M gal = 0.003, r nuc /r gal = 0.0002 and M bh /M gal = 

0.001. On a time scale of ~ f r (rbh) ~ 10 5 f ga ] the Bahcall- 
Wolf cusp is established near the black hole and the veloc- 
ity dispersion drops slightly due to the change in the density 
slope. At longer times the inward flux of heat from the galaxy 
begins to erase the dip in the velocity dispersion profile and 
causes the density to drop. 

In addition to heat transfer from the galaxy, destruction of 
stars by the black hole is effectively a heat source which by 
itself would cause the density of the nucleus to fall ([Shapiro 
119771) . However, integrations in which the loss term J (E,t) 
was set to zero evolved almost identically to those which in- 
cluded the loss term. Figure|7]shows an example. Apparently, 
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FIG. 9. — Evolution of the central density in the black hole models illus- 
trated in Fig. [8] The vertical axis is the density at a radius of 0.1 times the 
black hole's initial influence radius, compared with its value at t = 0. Left: 
time in units of galaxy half-mass relaxation time. Right: time in units of the 
galaxy's initial relaxation time evaluated at the black hole's influence radius. 

in all these models, the difference in temperature between nu- 
cleus and galaxy, and not the destruction of stars by the black 
hole, is the dominant effect which causes the nucleus to ex- 
pand. 

Figure [8] illustrates the early evolution of models from Se- 
ries I-III; the figure shows density profiles at t = 0, and af- 
ter an elapsed time of 10 4 f ga i. The top row of the figure 
shows the effects of varying the nuclear radius in galaxies 
with Mnuc/Mgai = 0.003 and fixed black hole mass (Series I). 
In each of these models, the initial velocity dispersion profile 
exhibits a dip outside s» 0.015, and the early evolution is 
a combined result of the inward flow of heat from the galaxy, 
and formation of the Bahcall-Wolf profile at r <, r\,h due to 
scattering of stars onto tightly-bound orbits. By ~ 3 x 10 5 f ga i 
almost all differences between the initial profiles have been 
erased and the nucleus has been transformed into ap~ r~ 7 / 4 
cusp with a mass approximately equal to the initial nuclear 
mass. 

The middle row of Figure [8] shows the result of varying the 
nuclear mass at fixed nuclear radius; the black hole mass is 
set to M nuc /3 (Series II). In low-mass nuclei, the evolution 
is similar to that of the models in Series I, with the differ- 
ence that the final normalization of the Bahcall-Wolf cusp is 
proportional to the initial nuclear mass. For M nuc ^ 0.01M ga i 
(Mbh <5 0.003M ga i), the initial velocity dispersion profile lacks 
a local minimum and the heat flow from the galaxy is greatly 
reduced; the evolution in these models is dominated by for- 



mation of the Bahcall-Wolf cusp. 

The bottom row of Figure[8]shows the effects of varying the 
black hole mass on a model with fixed (M nuc , r nuc ) (Series III). 
Here the differences are due almost entirely to the effect of 
the black hole on the velocity dispersion profile: larger black 
holes imply a smaller, initial dip in the temperature profile 
and less evolution. In these models, the contribution of the 
loss term to the changes in p(r) is most clearly seen, but it 
is still slight. It is interesting to note that the form of p(r) at 
r < rbf, can depart strongly from the "zero-flux," p ~ r~ 7 / 4 
form at early times in the models with larger black holes. 

Evolution of the central density over the full course of the 
integrations is shown in Figure [9] for the same models in 
Figure [8] Plotted there is the density at a (fixed) radius of 
0.1 rw,, normalized to its initial value. Aside from an early in- 
crease due to the formation of the ~ r~ 7 / 4 cusp, the general 
trend is for the density to decrease; the decrease is smallest 
(all else equal) for bigger black holes, for the reasons dis- 
cussed above. ForMbh/M ga i = 0.001, the density declines as 
p ~ f~ - 5 at late times (f > f,-(^/,)), with steeper(shallower) 
falloff for smaller(larger) My,- 

What do these models imply for the observed properties 
of nuclear star clusters? The general tendency of the black 
hole models is to evolve toward lower central density at fixed 
nuclear size, i.e. to move toward the left on Figure [T] and 
eventually to disappear from this plot as the nuclei become 
too diffuse to be detected. Evolution is strongest in the 
case of models containing nuclei of low initial density and 
small (Mbh SS 3 x 10~ 4 M ga i) black holes: these galaxies have 
a strong "temperature inversion," implying nuclear heating, 
while the presence of the black hole inhibits core collapse. 
The heavy curves in the lower left panels of Figures [8] and [9] 
show that the nuclei in such models are the most rapidly de- 
stroyed, and assuming the presence of such "small" black 
holes in nucleated galaxies would help to explain the lack of 
points in the upper left part of Figure [T] Larger black holes 
(Mbh i=J M nuc ) tend to stabilize nuclei by reducing the tem- 
perature inversion and decreasing the flow of heat from the 
galaxy. In such models, a Bahcall-Wolf cusp is established in 
a time ~ frt/bh) ~ 10~ 5 f ga i and little evolution subsequently 
takes place. 

While no attempt was made to construct evolutionary 
black-hole models that match the current state of NGC 205, 
it is clear from Figures [8] and [9] that including a black hole 
of mass Mbh <D 10~ 3 M ga i in this galaxy would cause the ini- 
tial configuration to evolve very little over 10 Gyr, aside from 
the establishment of the p ~ r~ 7 / 4 cusp inside ~ 0.2rbh ~ 
0.1 pc(Mbh/lO 5 M ). The kinematical upper limit to Mbh m 
NGC 205 is - 5 x 1O 4 M (IValluri et al J 120051) so the influ- 
ence of a black hole would be limited to ~ 0.05 pc w 0.02", 
below HST resolution. 

4.3. Results: Rates of Stellar Tidal Disruption 

An unambigu ous signature of a massive black hole is a tidal 
disruption flare (Ko mossa & B ade 1999). Figure [TOl shows 
the rate at which mass (in stars) is scattered into the black 
hole's disruption sphere for the various models. The left and 
middle columns of this figure show M expressed in units of 
M ga i/f ga i (which does not change with time) and Mbh/tr(rbh) 
respectively; the latter unit is itself a function of time via the 
time dependence of p(Vbh,?) and cr(n>h>0- (We emphasize 
that the mass of the black hole is assumed constant in these 
Fokker-Planck models.) The right column in FigureflOlshows 
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FIG. 10. — Left column: Rate at which mass is scattered into the tidal disruption sphere of the black hole, expressed in units of the galaxy mass divided by 
the galaxy half-mass relaxation time (both of which are constant). Middle column: Same as the left column, except that the loss rate is expressed in units of 
black hole mass divided by the (time-varying) relaxation time at the black hole influence radius fjj. Right column: The integral until time t of dM/dt, i.e. the 
total mass scattered into the black hole's disruption sphere, in units of the (constant) black hole mass. Lines and line widths correspond to the same models in 
Figures [8] and[9] The bottom panel shows, in addition, models from Series IV, with log 10 W = (6,7, 8,9, 10, 11); line width decreases with increasing N. 
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time [yr] 

FIG. 1 1. — Top: Stellar tidal disruption rate vs. time for models from Se- 
ries IV, which have fixed initial structural parameters but which vary in their 
assumed total mass M„ a i and hence in the number of stars, N = M„ a i/M ( :;. 
Curves correspond to N values ranging from 10 6 (thickest) to 10" (thinnest) 
increasing in steps of 10; the corresponding galaxy luminosities range from 
Mg — —8 to Mb — —20.5. Black hole masses are fixed at 10" 3 M gal . Dashed 
lines are defined in the text (eq. l36t . Bottom: Total disrupted mass. Galaxy 
relaxtion times have been assumed to scale with galaxy luminosity as in eq.[2] 

the total mass lost to disruptions, integrated from time zero 
until time t , in units of Mbh- 

While the variation of mass loss rate with time is com- 
plex, the general trend is for M to decrease as the nuclear 
density drops. We consider first a "fiducial" model defined 
by M nuc /M ga i w 0.003 and Mbh/M ga i ps 0.001; these are ap- 
proximately equal to the average mass ratios that are observed 
in galaxies with nuclear clusters and black hole s respectively 
(IWehner & Harrisll2006t iFerrarese et aT1l2006al) . We also fo- 
cus first on the case N = 10 10 which was assumed in Se- 
ries I-III; this corresponds to a galaxy absolute magnitude 
M B ~ -18 (eq. [25} and to r gal » 7 x 10 14 yr (eq. fl. By 
( = 2x 10~ 5 f ga i m 15 Gyr, the Bahcall-Wolf cusp has been 
established in these models and the loss rate to the black hole 
is essentially independent of the initial nuclear radius: 



M remains relatively constant thereafter if expressed in these 
units. The total mass in disrupted stars at t ps 2 x 10~ 5 f ga i in 
the same subset of models is ~ 0.02Mbh, implying a mean 
disruption rate over this time of ~ 10 3 Mbh/f ga i, or 

tin 1O- 6 M yr- 1 ( J**-) ( -^-Y 1 . (33) 
° 3 \ 1O 6 M J \ 10 14 yr J v ; 

However, much higher and much lower rates are also pos- 
sible: 

- At early times, before the Bahcall-Wolf cusp is fully estab- 
lished, M in the fiducial models is higher for smaller r nuc / r ga \ 
by as much as ~ an order of magnitude. 

- The second row of Figure[l0]shows that the initial M scales 
roughly with M nuc for fixed r nuc /r gal and M hh /M nuc . 

- The third row of Figure [10] shows that the mass loss rate at 
early times (i.e. before the nuclear density profile has been 
strongly affected by the presence of the black hole) is a de- 
creasing function of Mbh a t fixed (M ga i,f ga i) , for r easons that 
were discussed in detail by Wang & Merritt d2004l) . 

The results discussed so far have described the effects on M 
of changing the structural parameters of the initial galaxy. An 
equally important parameter is the assumed number of stars 
N. In addition to fixing the unit of time f ga i, N also deter- 
mines the details of the loss-cone refilling process and hence 
the number of stars scattered into the black hole in one relax- 
ation time. Small Af corresponds to a "full loss cone" and large 
N to an "empty loss cone." (Small values of Mbh also imply 
a full loss cone; for instance, in the third row of Fig. [TO] the 
full-loss-cone regime is reached forMbh/M ga ] ps 10~ 3 .) 

Recalling the definition of q(E) (eq. [27b as the dimension- 
less ratio between the orbital period and the loss cone refilling 
time, these two extremes correspond to q 3> 1 and q -C 1 re- 
spectively. It can be shown that the number of stars scattered 
into the loss cone in one relaxation time scales with q as 

7 xt«[lnj2/4] -1 , ? «1; (34) 

ocq- l ocN, #>1 (35) 
where the q °< N~ 1 dependence assumes a galaxy with fixed 
structural properties, i.e. fixed (M ga ],r ga i). This decrease in 
Mf ga i at small is visible in the last row of Figure [lOl for 
N <, 10 8 , i.e. M B k, - 13, the loss cones in these models are 
essentially full and M is independent of the relaxation time. 

The A^-dependence is illustrated more clearly in Figure [TT] 
which shows M and J Mdt for the same models, expressed 
this time in units of solar masses and years; equations (O 
and d25T > were used to convert Af and f ga i into physical units. 
Each of these models evolves in nearly the same way when 
expressed in terms of t/t sa \; as discussed above, changing the 
flux of stars into the loss cone has very little effect on the 
density evolution. The differences in the curves of Figure QT] 
are due primarily to the different time scaling (small galax- 
ies evolve faster) and secondarily to changes in q(E) as p(r) 
evolves. 

The initial disruption rate scales almost linearly with 
galaxy mass. This dependence is most easily understood 
in the small-AVfull-loss-cone limit, where M oc M ga i/? cr oc 
(M ga i/f r h) (frhAcr) with f cr the galaxy crossing time. From Fig- 
ure QJl, f r h ~ M ga i, and f r h/?cr ~N ~ M ga i, SO M — M ga ] — N. 

This is similar to the dependence rep orted in other studies, 
e.g. Figure 21 of iFreitag et al.l (120061) who considered black 
holes in power-law nuclei. 

After the Bahcall-Wolf cusp is established, the evolution of 
M is dominated by the systematic decrease in nuclear density. 
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The dashed lines in the top panel of Figure QT| which accu- 
rately reproduce the disruption rates for t J> f r ('"bh)j are 
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FIG. 12. — Points are stellar tidal disruption rates computed by Wang & Merritt (2004) for a sample of early-type, "power-law" galaxies, with luminosity 
profiles parametrized as "Nuker" laws. Open circles are galaxies for which the inner (projected) slope parameter T of the Nuker-law fit was 0.9 or greater, 
corr espo nding to a p ~ r -2 spatial density profile for the inner parts of the galaxy. Lines show disruption rates in the Fokker-Planck models from Series IV 
(Fig.ED at four times, 0.01 Gyr (thin line), 0.1 Gyr, 1 Gyr and 10 Gyr (thick line). 

might have been over-estimated, such a conclusion is proba- 
bly premature until detailed models can be constructed based 
on the newest luminosity profiles for these galaxies. Six of the 
galaxies in the Wang & Merritt sample are class ified as "nu- 
cleated" (Type la or Type lb) bv lCote et alj (|2()06). Inspection 
of the ACS luminosity profiles of these galaxies suggests that 
none of them exhibit a clear, two-component structure like 
that of NGC 205; a "Nuker" model might therefore provide 
an adequate fit to these galaxies. 

Figure [12] also emphasizes the strong time dependence of 
N in low-luminosity galaxies, due to the structural changes 
caused primarily by transfer of heat from the galaxy. This 
plot suggests that tidal flaring rates in dE galaxies are likely 
to be a strong function of the time since nuclear formation. 

5. SUMMARY 

1. Half-mass relaxation times in nuclear star clusters are 
interestingly short, ranging from ^10 Gyr in galaxies with 
absolute magnitude M B = — 18 to ~ 1 Gyr alM B = -15. 

2. In the absence of a massive black hole, evolution of nu- 
clear star clusters is a competition between core collapse and 
heating from the external galaxy. Observed nuclei appear to 
be in the core-collapse regime, with the exception of nuclei in 
the brightest galaxies whose relaxation times are so long that 
they would not evolve appreciably in 10 Gyr. 



The bottom panel shows time-integrated masses; these are 
roughly 10~ M bh at t = 10 Gyr in all the models. 

Tidal destruction rates computed in the past for lumi- 
nous E galaxies have generally been found to decrease 
with increasing galaxy/black hole mass (ISver & Ulmerl fT99l 
iMagorrian & Tremaine1ll999t IWang & Merritt! |2004|) 7 oppo- 
site to the trend shown here in Figure QT| One reason is 
probably the rather different models used in these studies for 
representing p(r), e.g. "Nuker" m odels which have unbro- 
ken power-law cusps at small radii (Bv un et al.| [l996). These 
are reasonable models for some galaxies but do not repro- 
duce the nuclear star clusters of low-luminosity ellipticals 
(Ferrarese et al. 2006b). Figure Q~2]plots disruption rates from 
Wang & Merritt] (120041) for a sample of galaxies modelled 
with Nuker profiles, compared with dN/dt from the Series 
IV models at four different times: 0.01, 0.1, 1 and 10 Gyr. 
The rates comp uted here at t = trace the lower envelope of 
the points from lWang & Merritt! (120041) . 

While these results suggest that t idal disruption rates for 
some of the galaxies modelled by Wang & Merrit j (120041) 



Nuclear Star Clusters 



13 



3. Adding a massive black hole "turns off" core collapse, 
and nuclear star clusters containing black holes always ex- 
pand, due primarily to heat input from the galaxy, but also 
due to the effective heating associated with stellar tidal disrup- 
tions. Nuclei with modest black holes dissipate more quickly 
than black-hole-free nuclei; large black holes, on the other 
hand, tend to slow the transfer of heat from the galaxy by re- 
ducing the temperature gradient. 

4. Rates of stellar tidal disruption generally decrease with 
time due to the drop in density of the nucleus. Rates are lower 
in galaxies of lower mass, assuming that the relation between 
black hole mass and bulge mass extends to low galaxy lumi- 
nosities. A simple expression is presented that reproduces the 



time dependence of the disruption rate in a restricted class of 
models. 
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namics. Stefanie Komossa also provided useful comments. 
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